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(57) Abstract 

A method and automatic system (Fig. 1) for extracting 
both ordinary and unique features of a stripe pattern like a fing- 
erprint from a gray level image wiihout binary processing. A pre- 
cise direction array (fig. 3) with th>: same number of image point? 
as the image is generated by calcinating the average direction of 
local ridges for every point in the i mage with a quick recurrent al- 
gorithm. There is also generated ? curvature array in which each 
element presents the accuracy of the local average direction at 
corresponding points in the direction array and image. A region 
of clear ridges extracted from bai Aground and noise before de- 
tecting image features. The riege trends (Figs. 8a-8e) and 
forkedness of a point axe decided yy analyzing the distribution of 
ridge directions on a circle around the point. Fox finding the 
cores and delta?, trend analysis is used for each singularity that is 
a maximum point on the curvanire array. The coordinate axis 
that is consistent For various fingc rprint types is decided by anal- 
yzing the structure of the direction array macroscopically. 
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IN GRAY LEVEL IMAGE 



10 



1. ftACKGRClUNH OF THE IN\ 

1.1 Field Of The Invention 
This invention relates to the automatic detection of both 

the common f eatures (i.e- cores, deltas and minutiae) and the 
unique features (shape and global features) o£ a fingerprint by 
processing a gray level image of the fingerprint* 

1.2 Description Of The Prior Art 
The history of identifying and verifying individuality 

according to dermatogliphic features is very long* Chinese 
people used to print their palms and fingers on documents and 
contracts as credit even since the seventh century A.D V 
Although there are many other methods for identifying 
individuals today, fingerprint identification is still the most 
15 widespread Jind credible. However, since fingerprints became a 
legal identifier of persons about one hundred years ago, the 
number of fingerprint records has grown very quickly and manual 
management of the files has become very difficult- As a result, 
many automatic and semi-automatic systems for processing, 
20 recognizing, searching, and identifying fingerprints have been 
proposed. 

The mo«t widely used method for detecting features in many 
present automatic fingerprint identification systems is based 
upon binary image processing. A binary image is one in which 
each image element has one of only two binary values, e.g. o or 
l. The key procedures in such processing are image enhancing, 
binarizing, thinning, smoothing and modifying. The minutiae of 
a fingerprint are detected by scanning the thinned binary image 
with a 3x3 window. Usually, the cores and deltas of finger- 
prints are also detected by scanning a binary image. A core 
exists where one or more lines of a fingerprint form a closed 
path, frequently a circle, or undergo an abrupt 180' direction 
change, A delta, as used herein, exists when three ridge lines 
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meet at a ciommon point. A delta may be more accurately 

referred, to as a Y. 

in U.S. Pat. No. 4,083,035, an apparatus is provided for 
detecting th* position (X and Y) and orientation angle (8) of 
5 minutiae in a binary data bit stream of a 256x256 thinned 
image. The mi-nutia orientation detector obtains an 8-bit vector 
average of aLl local angles present in each of a plurality of 
8x8 bit winehws across the image. This vector average of all of 
the local singles within a given 8x8 bit window is the 
10 orientation angle fl for each minutia that is positioned within 
that given 8:e8 bit window. There are 32x32 such windows on the 
Image, i.e. » 32x32 ridge orientation array will be generated. 

In U.S. Pat. No. 4,310,827, the minutia direction of an 
ending is defined as the direction of a single direction vector 
15 drawn from i:he ending to an arrival point, i.e. a skeleton 
point locate! by tracing a predetermined accurate length from 
the ending. The direction of a bifurcation is defined by a 
• direction symmetrical to the average vector direction of three 
arrival points. 

20 in U.S. Pat. No. 4,156,230, a 7x7 template scanning window 

is passed electronically over a 29x29 sub-array of the 32x32 
ridge contour data as in U.S. Pat- No- 4,083,035 to generate a 
set of correlation values corresponding to each contour data 
element and to a plurality of reference angle vectors. The 

25 correlation values are processed for determination of peaks and 
valleys. The resultant data, representing the number of 
correlation ;?eate and the direction of each, provides 32 values 
which define the location and angular orientation of cores and 
deltas of a fingerprint. 

30 in U.S. Pat. NO. 4,151,512, the topological data, 

identifying singularity points such as tri-radii (i.e. deltas) 
and cores, a* well as ridge flow line tracings related to those 
points, are extracted from a 32x32 ridge contour array as in 
U.S. Pat. No. 4,083,035. Subsequent to making the first cell 

35 tracing in any one direction from a tri-radius or core point, 
the information from the ridge contour array is used to supply 
additional angle data to continue each trace. Some logic 
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circuits determine the next row and column address incremental 
values according to a specification chart. The maximum length 
of a trace :.s 48 cells on the ridge contour array „ Based upon 
the number cf singularities located , an initial classification 
can he made wherein an arch is identified if no tri-radii are 
lgcated, a whorl may be identified if two tri-radii are located 
and a generM loop type may be identified if one tri-radius is 
located. Th« loop type pattern is classified according to the 
direction and size of the flow tracings by comparing them with 
a set of prc&stored references. 

The following publications are also of relevance to the 
present invention: 

Shen, "Several local properties of digital picture and 
their applications to the extraction of descriptive information 
of fingerprints" , Acta Scientiarum Naturalium r Universitatis 
PeJcinensis, NO. 3,. 1956, pp. 38-51. (in Chinese). 

Shen, "tfhe digital pseudo-curvature and its applications" , 
Applied Mathematics, Sept. 1988, No* 3, Vol. 3, pp. 382-391. 
(In. Chinese) . 

Shen et. al- r W A similarity Measurement and classification 
of Fingerprint", Proc. of 4th Chinese Conf. on Pattern 
Recognition and Machine intelligence, 1984. (In Chinese). 

1.3 Probleiis In The Prior Art 
25 The problems listed here below relate to the manner of 

identifying, or designating, cores, deltas and the shapes of 
fingerprints based upon ridge directions in the patents cited 
above. 

There are two problems in calculating ridge direction 
30 array: (1) Usage of the same value of ridge direction for every 
point in an 8x8 window will produce serious errors when the 
window is in a region where ridges curve significantly, (2) 
There is no measurement provided for representing the accuracy 
of each average ridge direction - 
3S In U.S. Pat. No. 4,310,827, the direction of a minutia 

depend on the arrival points. So the direction will be effected 
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if any arriv.il point can not he found or if the skeleton ridges 
are not smaath enough. 

There are four problems in analysing ridge trends.- (1) A 
7x7 window ii a ridge contour array, i.e. a Sfixse window in the 
5 original image, is too large to find cores and. deltas of snail 
whorls, or loops. (2) The fixed window size is not suitable for 
various types of cores and deltas. (3) The angular orientation 
with 32 values as well as cores and deltas with 29x29 positions 
are not accurate enough. (4) As many as 841 (-29x29) elements 

10 have to be analyzed for every fingerprint. 

There are three problems in ridge flow tracing: (l) Each 
step in ridje flow tracing passes 8 points because every 
element in the ridge contour array refers to an 8x8 region in 
the image. This is too large for tracing at regions where 

15 ridges curv«* significantly. (2) The errors of position and 
direction are not accumulated to correct the trace. (3) The 
next step mey be wrong when a core, delta or noise region is 
touched in tracing. 

There are three problems in classification: (1) The 

20 initial classification based upon the number of deltas may be 
wrong in case an existing delta can not be found. (2) The loop 
sub-olassifi cation by comparing the rough flow tracing is 
sensitive tc the initial fingerprint impression. (3) There is 
no sub-classification for whorls- 

25 Finally, the main factor which effects the accuracy of 

fingerprint features extracted by binary processing is that 
much original information in a gray level image of the 
fingerprint may be lost after binarizing. 

30 2. SUMMARY NP TTTR IHVEMTIOH 

2.1 Objects Of The invention 

It is therefore a. main object of the present invention to 
extract cores, deltas, minutiae, and shape and global features 
of a fingerprint from a gray level Image based on original 
35 information as much as possible. 

Xt is a more specific object of this invention to provide 
a. quick algorithm for calculating the average direction of 
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local ridges i at every point of a fingerprint and for generating 
a precise direction array* 

Anothe: : oh j ect of the invention is to provide a 
measurement . termed local curvature, representing the accuracy 
5 of each local direction that is easy to calculate. 

Another object of the invention is to provide a method for 
separating u region of clear ridges from background and noise 
in the imag<». 

Another object of the invention is to provide a method for 
10 analyzing the ridge flow trends around any point in the image 
to decide iks trend directions and forkedness. 

Another object of the invention is to find the cores and 
deltas of a fingerprint by analyzing the trends only for each 
singularity of the image rather than analyzing the trends for 
15 every point of a direction array* 

Another object of the invention is to .locate the position 
of the center and central orientation of a plain arch of a 
fingerprint, 

Another object of the invention is to establish a 
20 coordinate axis of any fingerprint that is consistent for 
various typos and shapes of fingerprints. 

Another object of the invention is to accurately trace 
shape lines, contour lines and normal lines of a fingerprint. 
Another object of the invention is to classify 
25 fingerprints! according to the structural relations among shape 
lines. 

Another object of the invention is to extract shape 
features from the shape lines that are consistent for both 
whorls and loops, and to further classify fingerprints 
30 according tc i the shape features. 

Another object of the invention is to extract global 
features of any fingerprint , including plain arch, however 
imperfect or partial it is and whatever type or shape it has. 
Another object of the invention is to calculate the global 
35 difference between two fingerprints to finely classify and 
•distinguish them. 
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Another object of the invention is to detect minutiae and 
their attributes from gray level images of fingerprints . 

Another object of the invention is to calculate both the 
quality level, and vector of fingerprints with regard to several 
5 aspects, for example noise level r area of clear region, 
position of center, number of minutiae, etc. 

2,2 BRIEF DlgS^ TWnflT OV fWR DRAWING 

FIG. l :ts a processing flow diagram showing the sequence 
10 of basic steps characterizing the present invention. 

PIG. 2 shows a point of an image and its neighborhood used 
to explain the calculation of local ridge direction and 
curvature according to the invention* 

FIG* 3 shows a point and its four adjacent points for 
15 calculating ;:our gradient models. 

FIGS- 4*-4e show a symmetric convex region and its four 
subsets for calculating four average gradient models. 

FIGS* 5i and 5b show the neighborhoods of two adjacent 
points and their common area as well as parts of the 
20 neighborhood!; of two adjacent points and their common points. 

FIGS- 6c.-6d show various octagonal regions of clear ridges 
in a fingerp:rint image • 

FIGS- Va-7d show various types of fingerprint core 

patterns. 

25 FIG- 7e shows a typical fingerprint delta pattern. 

FIGS, asi-Se show the ridge trends of various singularities 
and analysis circles. 

FIGS* *a-9e show the difference values within trend 

analysing, 

30 FIGS, ida-lOr show the shape lines of various fingerprints 

for 18 classes of shapes, 

FIG. 11 shows the macroscopic structure of peripheral 

ridges of a fingerprint. 

FIG- 12 shows a vault line and normal lines for locating 
35 the center o:? the coordinate axes of a fingerprint. 

FIG, i:i shows a manner of determining the central 
orientation of a plain arch. 
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FIG. 1A illustrates the extraction of shape features 
extracted on shape lines » 

FIG. lit illustrates the extraction . of global features of 
a fingerprint. 
5 FIGS* :L6a«16e show various minutiae; ' 

FIGS. I7a-17d show the basic features of minutiae in terras 
of a ridge or valley respectively near a core or delta. 

FIGS. 18a and 18b show two neighboring points in a 
tracing. 

10 FIGS * 19-24 are pictorial views illustrating various 

tracing operations according to the invention on gray level 
fingerprint images, where FIGS- 19, 20 and 21 depict the 
tracing of lines for locating the center and center orientation 
of fingerprint patterns containing a whorl, a loop and an arch, 

IS respectively; and FIGS. 22, 23 and 24 depict the detection of 
global features for an arch, a loop and a whorl, respectively , 
according tci the invention on the basis of the local directions 
of ridge lines of the fingerprint image patterns at points 
arranged along concentric circles. 

20 FIG* 2) is a pictorial view illustrating extraction of 

minutiae f rc<m the region of an arch in a gray level fingerprint 
image. 

2.3 GENERA1., DEFTKTTIOKS 

25 There are many constants, parameters, variables, and 

functions used herein. Some of these are defined in the C 
programming language as follows: 

x[] nutans an array named by x; it can be defined as a set 
in which 

30 x|]] - {x[0], . x[m-l]}, 

where m>0. 

x[] [1 means a two dimensional array, or matrix; it can be 
defined as a set in which 

x|;][] - {x[0]£G], x[m-l][n-l]) f 
35 where m>0, n>0. When each member of such a matrix is of the 
type x[i][j; , then, for each member, i is called the line, or 
row, number and j is called the column number. 
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x = C ? y : z; means that if condition C is true, t hen 
x^=y? else as- 
sign (x) " (x <o)? -1 : l; 

int(x) means the largest integer that is not greater than 
5 Therefore, if x>0, then int(xH) ♦ 5) means the most 

approximate integer of x. 

x%y mefLne the remainder, that is not smaller than 0 f when 
x is divided by y. 

(x,y) :aeans a digital point at coordinates x and y of a 
10 plane. It cim also mean, in the appropriate context, a vector 
from origin ;?oint (0,0) to (x,y> . 

dv(X) maans the direction of a vector X, which direction 

is in the raige [O, 

#A means the cardinal or the number of members in set A. 
23 S(y(X) f X, A) means the sum of values y(X) for every 

member X in iset A. 

|xj means the model or absolute value of x. if x is a 

number then 

|x[ => (x<0)? -x s x; 
20 if x^Cx,,, ..,x n ) is a vector, then 

|x| = ^EUl,i,{l,...,fltt; 
p_x mee ns a parameter which may be predetermined or 
calculated for use in performing computations according to the 
invention. The range of each predetermined parameter will may 
25 be used in preferred embodiments of the invention will be 
listed below, 

w = 3,1<1153; 
p_jr = 2!J2. 

30 p pBCTTPTlOM OF TP * P^FE RRED EMBODIMENTS. 

3.1 General Procedure 

Referring to FIG.l, the general procedure of this 
invention foe extracting as many features as possible of a 
fingerprint :!rom a gray level image without binary processing 
35 is shown. 
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The significant steps in an exemplary method according to 
the invention are: 
* 1) Input and Pre-processing 

The in; put digital image has L rows and K columns of image 
5 elements, els shown in FIG.l. The intensity , t or brightness 
level, in each image element has 3 to 1024, and preferably 256, 
gray levels and the image element density is 500 dpi (dots, or 
image elements f eer inch) in each coordinate direction. For the 
further processing, the original range of gray levels of the 
10 fingerprint image is transformed into a uniform range ♦ 



2) Calculating Direction Array and curvature Array 

A dire<rtion array and a curvature array are calculated for 
15 each image )>oint by a quick recurrent algorithm, possibly with 
the aid of some tables. The direction array and curvature array 
values for an image represent the average ridge direction and 
its accuracy or variance at each point of the image* 

3) . Cleaning Background 

20 segmemijig an equiangular, or regular, octagonal region of 

clear ridges of a fingerprint from the background by eight 
straight lines according to the curvature array, and excavating 
all large conjunct regions of noise points in the octagon* All 
further processing is within this region. 

25 4) Trend Analyzing 

The ridge trends and forkedness of any point can be 
determined by analyzing the distribution of ridge directions 
around it o:z circles with different radii. 
5) Finding Cores and Deltas 

30 All singularities are found by scanning the curvature 

array with .* 3x3 window to find all maximum curvature points. 
All cores and deltas are located by analyzing the trends for 
each singularity. 

• 6) Line Tracing 

35 Contour lines, shape lines and normal lines of a 

fingerprint are traced , based on the direction array , 
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accurately by accumulating the errors of coordinates and 
directions t>i correct the trace point by point. 

7) Deciding the coordinate Axis 

By macroscopically analyzing the structure of a direction 
5 array, the fingerprint is centered on a coordinate axis system* 
meanwhile the central orientation of the fingerprint is 
selected on she basis of the trend at the center. 

8) GLissif ication 

A fingerprint is classified into one of 18 classes 
10 according to the structural relations among the shape lines. 

9) Extracting shape and Global Features 

Por classifying and fast searching, a few shape features 
are extracted consistently for various fingerprint classes 
except plain arch from shape lines. Furthermore a plurality of 
15 global featares are extracted consistently for various 
fingerprint classes with reference to the coordinate axes. 

10) Defecting Minutiae 

Based on the description of minutiae when the shape lines 
are valleys, all minutiae are detected by tracing each valley 
20 of the fingerprint in the gray level image. Each minutia is 
represented Isy its x, y coordinates and direction 6. 

11) Quality Checking 

A synthetic quality level and quality vector of 
fingerprints is presented based on the position of the center, 
25 number of minutiae, noise level, area of dear region, etc, in 
order to decide automatically, or to suggest to the operator, 
whether to accept, reject or reevaluate the fingerprint. 



30 



3.2 Direction Array and curvature Array 

All features of the fingerprint, including ordinary 
features (coi.es, deltas, minutiae) or novel features (shape and 
global features) are referenced to the average direction of 
ridges in a small region of the fingerprint. The calculation of 
local direction is very Important for image processing and 
35 feature extraction of fingerprints, in this invention, an array 
whose every element represents an average direction of the 
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textures in a small region of the image is. called a direction 
array * 

* There *re three features of the method provided herein for 

calculating a precise direction array; 
> 5 F;.T3t, the direction array is calculated directly 

from the gr*iy level image of the fingerprint, so the original 

information will be used as much as possible. 

Second r the direction array is calculated point by 

point in the image, i.e. every element in the array is a local 
10 average dir action of just one point in the image that is 

calculated on the basis of a neighborhood of the point, this is 

necessary spatially for regions on the fingerprint where the 

directions cf ridges change greatly or the curvatures are very 

high. 

15 Third, to represent the accuracy of the local 

direction al: each point, a local average curvature of the 
textures in the neighborhood of the same point is also 
calculated, and a curvature array of the image is generated 
therefrom* Tie usage of local direction at a point should refer 
20 to -the local curvature as the relation between an average value 
and its variance. A local curvature in some regions of a 
fingerprint, for example at a core, delta, scar or noise, will 
have very high curvature values which signify that local 
direction is meaningless there due to inconsistency of ridge 
25 directions and/or indistinct textures in that region. 
Generally, lawer curvature values mean better accuracy of the 
direction value at a point. 

For a civen image area S, which may include the entire 
fingerprint image or any selected portion thereof, there are 
30 derived four summation gradient values v f (i=l, 2, 3, 4), each 
representing the sum of the absolute values of the difference, 
g f (X), in gray scale image point values, f (X) , between each 
pair of points, X and X-Q lr in S for which Q, is a selected 
vector. Thus, as will be seen from the example to be described, 
35 depending on the value of v, is representative of the degree 
of change in the image across the image area in direction Q { . 
Moreover; v r will be particularly relevant to the point P at 
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10 



30 



tHe center o:r S because there is the greatest probability that 
it is the gradient condition at P which is described by v s . 

Table 1, below, provides an exemplary gray scale value 
matrix representing the gray values f (X) at respective points, 
X, of an image area S. Here, each point X has a horizontal 
coordinate k and a vertical coordinate 1. As is apparent from 
Table 1, the origin of the coordinate system will be somewhere 
above and to the left of the illustrated image area. 

Each po:Lnt X is represented by a pair of coordinates k, 1. 
In Table 1, the k and 1 coordinates at the center of s are n 
and m. 

Table 1 

s k| n-5 :a-4 n-3 n-2 n-i n n+l n+2 n+3 n+4 n+5 





m-5 {13 L4 15 15 IS 14 12 11 9 7 6 

m-4 J 14 L4 14 13 12 10 6 4 2 2 3 

m— 3 { 17 LI & 6 6 4 3 3. 4 6 8 

m-2 {10 7 4 



20 m-1 | 5 4 



3 



122 3469 11 
0 2 4 6 8 11 13 14 



3 4 5 6 8 9 11 12 14 15 15 
n+l j A 6- 9 10 12 14 16 15 IS 14 12 
m+2 ! 7 9 12 13 14 14 IS 14 13 10 8 



m+3 { 13 14 14 14 
25 m+4 I 15 L5 14 



14 13 10 8 8 6 4 



m+S | 14 12 10 



13 11 8 3 2 5 5 5 
87401578 



For each vector Q { there is produced a set S,- containing all 
points X foi- calculating g.(X) = | f (X) -f (X-Q«) | . The points in 
each set S, are obtained by using the corresponding values for 
Q |r i.e., s, is obtained by using Q,, S 2 by using Q a , ate. The 
values for g, (X) in each set S, are given in Tables 2-5, below, 
where hav« the following k,l coordinate values, 
Q t = (1,0)? 
35 Q 2 = (1j 1) ' 

Q, = (0,1)? 
=(-:.,!). 
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Table 2 





s 1 kj n-5 
1 | 


n-4 


n-3 


n-2 


n-1 


n 


n+1 


IH-2 


n+3 


TH-4 


n+5 




m-5 ! 


l 


1 


0 


0 


l 


2 


1 


2 


2 


1 


5- 


ln-4 } 


0 


0 


1 


1 


2 


4 


2 


2 


0 


1 




m-3 | 


2 


2 


3 


0 


2 


1 


0 


1 


2 


2 




* m-2 ! 


3 


3 


3 


1 


0 


1 . 


1 


2 


3 


2 




m-l { 


1 


1 


3 


2 


2 


2 


2 


3 


2 


1 




* [ 


1 


1 


1 


2 


1 


2- 


1. 


2 


1 


0 


10 


mfl | 


2 


3 


1 


2 


2 


1 


0 


0 


1 


2 




m+2 { 


2 


3 


1 


1 


0 


1 


1 


1 


3 


2 




m+3 j 


1 


0 


0 


0 


1 


3 


2 


0 


2 


2 




m+4 | 


0 


1 


1 


2 


3 


5 


1 


3 


0 


0 




nH-5 j 


2 


2 


2 


1 


3 


4 


1 


4 


2 


1 



15 



20 



25 



30 



Sj k 
„1_ 
m-5 
m-4 
m-3 
m-2 
m-l 
m 

m+1 
m+2* 
. m+3 
m+4 
m+5 



n-5 n-4 



Table 3 
iv 3 n-2 n-1 n 



n+1 n+2 n+3 n+4 n+5 



1 


0 


2 


3 


5 


8 


B 


9 


7 


4 


3 


5 


B 


7 


8 


7 


3 


0 


4 


6 


6 


7 


8 


4 


4 


1 


1 


3 


5 


5 


6 


4 


4 


1 


2 


4 


5 


7 


7 


5 


1 


1 


3 


8 


7 


7 


6 


6 


4 


2 


3 


5 


5 


6 


6 


6 


4 


3 


0 


3 


5 


6 


4 


4 


2 


1 


1 


2 


5 


6 


7 


5 


2 


1 


1 


4 


7 


6 


7 


6 


2 


0 


1 


3 


6 


10 


8 


3 


3 


1 


3 


5 


6 


6 


7 


8 


2 


3 


2 


3 
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Table 4 



10 



IS 



20 



25 



30 



35 



S3 K l 




j 1 ** 


n— 3 


n-2 


n-1 


n 


n+1 


n+2 


n+3 


n+4 


n+5 


_± 1 
























m-5 { 
























m-4 } 


1 


0 


1 


2 


3 


4 


6 


7 


, 7 


5 


3 


m— 3 | 


X 


*> 


— > 


7 


6 


€ 


3 


1 


2 


4 


. 5 


m-2 [ 


3 


4 


5 


5 


4 


2 


0 


1 


2 


3 


3 


m-l f 


5 


3 


1 


1 


0 


2 




^« 


K 
■J 


4 


3 


m j 




0 


2 


6 


6 


5 


5 


4 


3 


2 


1 


JDnrx j 




2 


4 


4 


4 


5 


4 


3 


1 


1 


3 


UttjI. j 




3 


3 


3 


2 


0 


0 


1 


2 


4 


4 


. to | 
BH-J | 


D 


K 


2 


1 


0 


1 


5 


6 


5 


4 


4 


— « * 1 






0 


1 


3 


5. 


7 


6 


3 


1 


1 




-L 




4 


5 


4 


4 


3 


1 


0 


2 


3 












Table 


5 














n ^ 


il ** 


XI— 3 


n-2 


n-l 


XI 


n+1 


n+2 


n+3 


n+4 


n+5 


1 1 
























J- j 
rn — 5 ! 
























m-4 { 


0 


1 


1 


2 


2 


2 


5 


5 


5 • 


4 




m-3 { 


1 


3 


4 


6 


4 


2 


1 


1 


2 


3 




m-2 { 


1 


2 


2 


5 


2 


1 


0 


0 


0 


1 




m-l | 


2 


0 


2 


2 


0 


1 


2 


2 


2 


2 




m { 


1 


1 


5 


4 


4 


3 


3 


1 


1 


1 




m+1 } 


0 


1 


3 


2 


3 


3 


3 


l 


0 


1 




m+2 j 


1 


0 


2 


1 


0 


1 


0 


1 


1 


2 




m+3 } 


4 


2 


1 


0 


0 


2 


4 


5 


2 


2 




m+4 { 


1 


1 


0 


1 


2 


2 


5 


6 


1 


1 




nH-5 | 


1 


2 


3 


3 


1 


1 


2 


4 


0 


2 





The nimiser of values for g { (X) on each set S, is less than 
the number oC points in area S because each value for g f (X) is 
calculated oily for a pair of values of f (X) and f (x-Q,) which 
are both in i»rea S. Thus, for example, on set S, there will be 
no value of g,(X) for which the k coordinate of P is n-5 
because the Jc coordinate of X-Q 1 is n-6, which is outside of 
area S. 
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Then, for each set S f/ there is derived a value v ? where 
v f = sg^P) for ail values for g^P) on set S f . Where i = 1, 2, 
3, 4, the VcLLues for v f in the example shown in Tables 1-5 will 
be ^^167? Vg=437; Vj=337 and V 4 =194. 
5 From these values for v { , there will be derived four 

further values u ir as follows: 

u, = sulx (v v v 3 ) 
^ = nuix <v 2 , v 4 ) ; 
10 u 3 - mi .n (v, , v 3 ) ; 

u 4 = min (v 2 , v 4 ) . 

In the case of the example in Tables l*-5, ^=337; u 2 =437; 
Uj=167; u^lM. 



15 



These values are then used to derive: 



20 For thit example shown in Tables 1-5, e^36, 148865. 

The derived values for v, u and e. are then used to derive 
a local direction value, d(f , P, B) , and a local curvature 
value, c(f, P, S), 

Referr.\ng to Figures 3, 3 and 4 of the accompanying 
25 DRAWING, according to the 1986 and 1988 publications of Shen, 
the formulas for calculating the local direction d and 
curvature c at each point P at the center of a set S in image 
f are as follows: 

d(f,P,St) = sign (v 4 -v 2 ) • arctn ( (v r e)/(v 3 -e) ) . (l) 

1$!^^ <2) 
(Uj+u,) •<u z +u 4 )-2u s u 4 

30 
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Xn the example shown in Tables 1-5 , which example is 
associated with P(n,m), d = -23. 5* and c - 0.0927. 

The noraal values of direction d are limited between -f/2 
(upward in FCG- 2) and +1T/2 (downward in FIG. 2) . Here d=0 
5 represents the horizontal direction which points toward the 
right Table l and Fie. 2. According to Equation (2), the 
normal values of curvature c are all between 0 and 1. 
Specially/ c=0 represents the texture with a plain local 
curvature, c==l represents the texture with an abrupt local 
10 curvature. Additionally, for each background point or a noisy 
point P, botli of c and d will be set to a special value that 
c(P)«*255 and d(P)=255. 

S(P) is a neighborhood of P, it is convex and symmetric 
with respect to P and with respect to the directions 0* and 
U5 45" , respectively. For example, a digital square, disc and 
octagon with. P as center are all neighborhoods of this kind. 
Each S f is a subset of S by deleting some border points as 
shown in FIG-4 where '.'is in S and S f , •*• is not in S, 
(1=1,2,3,4). 

20 Nov ref.irring to figs. 5 , because there is tremendous 

complexity fcr calculating the direction array and curvature 
array per point, a quick recurrent algorithm with some tables 
is proposed ta reduce the complexity on the basis of three key 
points: 

25 First of all, because each gradient model |f (P)-f (P-Qf) [ 

will be used as many times as the number of the points whose 
neighborhood includes both P and P-Q f , four arrays o£ gradient 
models at four directions as in fig. 3 are calculated firstly 
and signed by g ; respectively: 

30 g,.(P) = lf(P) ~ f(P-Q>)|r where i=l, 2, 3, 4. (3) 

second, there are many common points in the neighborhoods 
of both P and its adjacent point P". So v.(P') can be 
calculated recurrently from v ( (P) by subtracting g,(X) for each 
X at the left side L(S f (P) of s f (P) and adding g,(X) for each X 

35 at the right iside R(S,(P)) of S f (P") , as shown in FIG. 5a where 
each is in S(P), each is in S(P") and each -O' is in 
both sets: 
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V f (f,P ,S S (P-)) - ^(fpFrSfCP)) - E^W^s^P))) 

+ SCg^Xj^X^RtS^P'))). (4) 
For example, after v has been calculated for P (n,m) , v i 
may be calculated for P (n+1, m) by subtracting from P (n,a) 
5 the values of g, (n-4, 1) , and adding the values of g { (n+6, 1) p 
in each case 1 taking on each value from m-5 to HH-5. 

Furthermore as shown in FIG, 5b , (where ■*■ is in L, 
is in L A ard "0« is in both of them, S(g,(X),L*) also can be 
calculated recurrently: 
uo S(g f (x; ,x,L-) - z(g,(X),x,L) - g,(P') + g i (P"). (5) 

Where ?■ is the top point in L and P n is the bottom point 
in L A shown in FIG« 5b. 

Third, some look-up tables are used instead of a series of 
calculation 3 for direction d and curvature c* 
15 To calculate the value of d, a table 

Ti - {td[0], td[l], td[Md]}; 
is created, where each term Td[k] in Td is defined as 

M[k] = arctan(VMd) f for k = 0, l r - /Md. (6) 
where Md is a predetermined integer. So the values of terms in 
20 Td .are all between 0 and tr/4. For any group of v f and e with 
the condition v 1 £ v 3 , let 

di£f_d=Td[:Lnt(Md- ( (v^e)/ (v 3 -e) J+0.5) ] - aretan( (v^e)/ (v 3 -e) ) . 
Then according to the continuity of the function arctanf ) , the 
value |diff_d| will be very small if Kd is large enough. So 
25 that by the table Td, when v^Vj, Equation (1) can be 
transformed approximately tos 

d(P) = sign(v 4 -v z ) - ((v t <v 3 )? 

TdrintCMd-CCv^e) ■ (v 3 -e) )+o,5] i 
(JT/2-Td[int(Md- ((V 3 -e) ■ (v r e) J+0.5) ]) . (!') 
30 The integer Md can be selected large enough to assure 
sufficient .accuracy for y. 

To calculate the value of c, another table 

T<3 - <TC[0J[O], TC[0][1], . .., Tc[MC][MC]}? 
is created vhere each term Tc[i][jJ in Tc is defined as 
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where Mc is a. predetermined integer. So the values of terms in 
To are all between 0 and 1. ror any group of Uj, let 

5 diffc - TcElnt(tto-u 3 / u l +0 * 5 )Hint(Ifc-u < /u 2 +0.5)] - 

(1^+113) * (1^+1^) -2 *u»*u 4 

Then according to the continuity of the function in Equation 
(2), the value |diff_c| will be very small if Mfc is large 
enough. Sc that by the table To, Equation (2) can be 
XO transformed upproxiJnately to: 

c (P) = To [int(«e*u3/Uj.+o . 5) ] tint Ud&^/Ug+o .5)1. 



The integer He can be selected large enough to assure 

sufficient accuracy of the formula. 

The lool^up tables may be provided, e.g. in a non-volatile 
15 addressable memory, where each entry corresponds to a 

respective value of Td or Tc. 

If the neighborhood for calculating the local direction 

and curvature is a square window with 2*r+l points as both its 

length and w:.dth, let L and K represent the numbers of rows and 
20 columns, respectively, of an image array, and the algorithm is 

as follows: 



<1> calculate g { 

g,[ij -| f ti] m-mi 1 f 

25 lo-=l,...,K-l; 

g 2 [l] CW- 1 *CX3 Ek]-«[l-i] C*- 1 ! I > 
k-l r - - - ,K-1? 

k=o, . . . 

30 g 4 [lKk:|-=|f[lKlC]-fCl^Hlc+l]|, for 1-1,. ,Ir-l, 

]c=0 , — , K-2 ; 



for 1=0,.'.. ,£-1, 

for l*=l, • ■ . 

for 1=1, . . - i l-1 r 
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goto <2>; 

<2> Accumulating g. of each column k on a vertical line 
{(0,k), ... (2-r,k)}, 

w^kj^tg^l]^]), 1, {0,l,...,2'r) r for k=*l, . - „ ,K-i; 
5 v 2 [k]=3:(g 2 [l][kj), 1, (l f 2...,2-r)) r for k=i M ..,K-x; 

nfc[kj=a:(g 3 [l][k]), 1, {i,2,...,2-r}), for k=l, . . . ,K-1; 

w 4 tk]-3:(g 4 [l][k]) . 1/ { (1,2.,, ,.2-r}), for k=0,. ,K-2; 
where k is 'the number of a column, l is the number of a row 
and varies from 0 to 2-r, each w f -[k] is a summation of g f for 
10 l^o,,,., 2 »r and assigned by k* 

goto <3>; 

<3> Accumulate (i=»l, 2, 3, 4) in current window in each 
g f to oh t air i v r 

v^sfw^j]), j, (l,...,2-r)); 
15 v 2 =2(w 2 [jj), j, {l,.„,2»r})r 

v 3 =S(w 3 [j]), j, (0,...,2-r)); 

v 4 =2(w 4 [j]), j, {0,...,2.r-l}) y 

goto <*l>; 

<4> Calculate curvature and direction of every point (k, 1) 
20 in current window and shifting the window to right, 
if (v T +\ 2 +v 3 +v 4 <p_v) f then 
c[l]rk|=255; d[l][k]=255; 

else c^JCkl^ctintCMc-u^^+.SJirinttMc-u^/Uj+.S)]; 
i:f ( v 3 ««v 1 ) , then d [ 1 ] [k]=sign (v 4 -v 2 ) • w/A ? 
25 e!Lse (d[l}[k] = sign(v 4 -v 2 ) - ( (v^v,) ? 

Td[Md- (v,-e)/(v 3 -e)+.5]: 
(ar/2-Td{Md- ) v 3 -3) * (v^eJ+.B] ) ) ; 

> 

) 

30 -k=k+i; 

if (k>K-r) goto <5>; 
- v^+w, [k+rj-w, [k-r] ? 

v 2 =v z +w ;! [k+r]-w 2 [k-r] ? 

v 3 =v 3 +w :; f)c+r]-v 3 [k-r-l] ; 
35 v 4 =v 4 -w, t [k+r-l]-wjk-r-l] ; 

goto <i\>i 
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<5> Put the current window at the beginning of the next line 
and recalculating w t : 
1=1+1; 

if tl£L-r) then return? 
5 w^k^tlifcj+gt [14X-] [kj-g^l-r-1] [Jc] , for k~l, . . . ,K-1; 

v z Ck]=w 2 i:k]+g 2 [l+r] [k]-g 2 [i-r] [k] , for te»l, . . . ,K-1; 
v 3 [k3=w 3 i:k]+g3[l-*^] Ck]-g 3 n-r] [k] f for k=0, . - . ,K-1; 
v 4 [k3^w 4 |.k]-tg 4 [l+rjtk]-g 4 [l-r3 [k}, for fc=0, . . . ,K-2; 
p_v is a threshold which is directly proportional to r 2 , 
10 and the condition v 1 +v 2 +v 3 +v 4 <p_v at a point means that the 
contrast of gray level in its neighborhood is very low- 
Usually direction at those points are ignored,, and the 
curvatures axe assigned a special value of 255. 

Here each term w^k] is associated with a column number 
15 k- For example, referring again to the values in Tables 1 
and 2 above, let r=5, then, 

w,Eu]Hi:(gl[lUn], 1, {m-5,m-4 r . . • ,w+5 J )=17/ 
similarly, ff, ;n-hl]=26. 

Now, an example is provided in Table 6, which is derived 
20 from Table 4, for r=2. The values of w 3 [k} (k=*n-5,n- 

4 r ...,n+5) are calculated as in Table 6, while Table 7 shows 
the values of v 3 at the points from line m-3 to line m+3 in 
Table 4* 



PAGE 39/82 ' RCVD AT 6/29/2005 1:40:37 PM [Eastern Daylight Time] 1 SVR:USPT0-EFXRF-1/3 1 DNIS:8729306 ' CSID:631 549 0404 < DURATION (»$$):2W8 



06/29/2005 13:53 FAX 631 549 0404 



STRIKER & STRIKER 



-> US PTO 



© 040 



WO 93/07584 



PCT/US92/OS446 



21 



w 3 k'| 



Table 6 

n-fi n-4 n-3 n-2 n-1 



n 



n+1 n.+2 n+3 n+4 n+5 



10 



m-3 j 
m-2 J 
n-1 j 
m j 
m+l J 
hh-2 j 



10 10 12 15 13 14 12 13 16 16 14 

l.L 10 13 19 16 15 11 10 12 13 12 

II 9 12 16 14 14 12 12 11 10 10 
llL 8 10 14 12 12 12 12 11 11 11 
II! 10 11 14 12 11 14 14 11 11 12 

III 11 9 9 9 11 16- 16 11 10 12 
1!2 12 9 10 9 10 15 14 10 11 12 



15 



For example.* 

w 3 [m+2 J [n]=tt 3 [i[H-l] [n]-g 3 [m-2] [n]+g3[m+2J [n] =14-2+0=12. 

Table 7 



v 3 kj 

1 j 


n-3 


n-2 


n-1 


n 


n+1 


n+2 


iH-3 


m-3 J 


60 


64 


66 


67 


68 


71 


71 


m-2 { 


6S» 


73 


74 


71 


64 


61 


58 


} 


6;: 


65 


66 


68 


63 


59 


55 


ID | 


5£i 


56 


60 


62 


59 


58 


57 


BH-1 J 


55' 


53 


62 


65 


62 


61 


62 




50 


49 


54 


61 


63 


64 


65 


m+3 ! 


5£: 


50 


53 


58 


58 


60 


62 



For example , 

v 3 [m] [n-l]«v 5 [m] [n-2]-w 3 [n] [n-4]+w 3 [m] [n+1] =73-10+11=74 . 
30 So ths.t with the above algorithm, only 4 additions and 

subtractions, are needed on average for calculating v 3 at each 
point. 

3.3 Removing Background 
35 For most fingerprint images, there always are noisy 

textures and. other features in background * Before extracting 
the salient features of the fingerprint, it is necessary to 
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segment a region of clear ridges, or valleys, from background 
and noise, i. e. to decide the position and boundary of the 
region of ol<iar ridges- The curvature value of a point in 
background iis always very high due to the low contrast or 
5 noise, so the clear region can be obtained by cutting off the 
points with high curvature values. 

Bef erring to FIGS* 6, the method proposed here is for 
locating an cigui angular octagonal region of clear ridges of 
the fingerprint. The eight edges of the octagon are all 

10 straight line is with predetermined equispaced angular 

orientations, e.g. CT , 45% SO* , 135', 180*, 225% 270" and 
315*, respectively. 

Bef erring to Fig* 6a, the edges of an equiangular octagon 
are assigned E u E a respectively, and the comers are 

15 assigned by the coordinates (x^y,),..., respectively - 

So each edge can be represented by an equation of straight 
line as folltws, 

or x+-y=*2+y 2 * 
or x-y=x 4 -y 4 ; 
or x+y=x^+y 6 ? 
or x-y^xa-y a ; 



equations, i.e. by eight parameters y 2 , x 4 , Yt* *6i Y&* 

x s , y B }. The octagon is obtained by cutting the image with 
eight lines according to the curvature array sequentially and 

30 described by only a bytes of the positions of eight lines. 
There are several shapes of octagons shown in FIGS. 6. The 
algorithm to locate eight edges is as follows: 
<1> Set the original edges of clear region such that 
j^=r; x^L-r-1? y 6 -r; y 2 =K-r-l; 

35 calculate the average curvature of every row and column 

in the currext clear region, and store them in the two arrays 
ad [ ] and ac*.[ ] respectively , i . © - 













Y=1V 


20 




























x-i-x 7 -y 7 ; 


25 
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acl[y]^S(c[y] [x] ,x, {x 4 ,x 6 +l, * ♦ . ,Xg)) , 

for y*=y«, y 6 +l, . y 2 ; 

acl[x]?=x(c[y] £x] ,y, {y 6 ,y 6 +i, . . • ,y 2 }) , 

for x==x 4 , x^+l, .. w x fi . 

5 lei: n=0? 

<2> if (n>p_n t ) then goto <5>? 

where pji, :Ls a predetermined l imitation , i.e. the times for 
cutting the image is not more than p_n 1 «. 

I*et Xq being the number or coordinate of .the row of the 
10 current array with the minimum average .curvature value, i*e., 
acl[y 0 ;=rain(acl[y 6 ] , ... r acl[y 2 ]); 
similarly lat y c satisfy , 

ackCx^^inCack^], . ack[x a ]). 
if(ack[x 0 ]/(x B -^+l) £ acl[y 0 ]/(y 2 -y 6 +l)) then 
15 goto <3>; 

else guto <4>; 

<3> Cut th<» area by horizontal lines, i.e. determine edges 
E z and E 6 . Let 

P-c_y^>_c_p'S(acl[yj, y, {y 6 , . . . *y 2 })/ ty z -y 6 +x) 
20 l=max{y | (acl[yj > p_c_y) & (y*y 0 ) & (y*y 6 ) } .; 

l z =*min{y | (acl[y] > p_c_y) & (y£y D ) & <y<y 2 ) ) ? 
i.e. 1 is the largest row number satisfying l<y 0 , l>y d and 
acl tl]>p_c_jr; and 1 2 is the smallest row number satisfying 
l>y 0 , lj<y 2 end acl[l]>p_c_jy. Then recalculate y 2 , y 6 and 
25 ack[} in the; current clear region as follows, 

ack[x]-z:<c[y][x], y f {y & , y 6 +l, . ,y 2 )> ; 

for x = x^ f x^+l.| • . * § Xg * 

let n=rt+l; 
30 goto <S>; 

where p. c. p is a predetermined parameter for calculating the 
parameters p_c_x and p_c_y. 

<4> Cut area by vertical lines, i.e. determine edges B 4 and 
E ft . Let 

35 P^c,jc=E_c_p'E(ack[i], i, {x^ . . . f Xg})/ (x^-x^+l) ; 

kf*xtax{'X | ack[x]>p_c_x & x^Xq & x£xj ; 
k 2 ^min(:t | ack[x]>p_c_x & & y*Xg>; 
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i.e. k, is tie largest row number satisfying and 
aoX[k,]>P.-C_*; and is the smallest row number satisfying 
1^*0, k^xg aad ack[k2]>p_c_x. Then recalculate 3%, x* and 
aekt] in the current clear region as follows, 

5 

acl[y]^53tc[yj[x] r x, {x^+l, - • - »Xa>> > 

for y=y 6 , y s +i, • - • ^ 

let n=n-i'l; 
goto <2>t 

10 <5> Cttt are£L by hypotenuse lines, i-e. determine the edges 
E 1( . Ej, E 5 anc. E 7 of the octagon enclosing the clear region or 
determine thn numbers Xj, Y 4 , a»d y 8 . Let 

P_c-_z^(p_c_x + p_P_y)/2? 
then the four numbers can be calculated as following, 
15 x^intz | Cao3[s]>P-C_ a ) * & (z>(x 4 +X a )/2) > ; 

where for z=i'^, Xg-1, (x^x^/^, 

ac:i[z] =S(c[l][k], (k,l), AgCz))? 
where (2) = Uy 2 ,«K (]&-!# *«.)# (Y2-XB+ Z r x a> } . 

a% | (ac3[z]>p_c_z) & (z^ 4 ) fi 

20 where for 8=:^, x 4 +l» •••/ (x 4 +Xg)/2, 

ao;>E8]=£Cc[lHJc], Ck,l>/ A 6 (2>)' 
where ^(z) ■• { <Y 6 r z) * (Y*«r«-1> r • • • ' (y 6 -V z » x <> > • 

y A = mir{a I (ac3(z]>p^c_z) & (z£y 2 > & (z^+Ya)/ 2 ) > ' 
where for z^r v 2 2 -l, • (Yz+Yg}/ 2 * 
25 ac:j[z]=S(c[l][lc], (k,l), M*)> ; 

where (z) - { Cy 2 , z> , (y E -i/ Z" 1 ) ' • • • ' Cy fi +^-z r * A ) > . 

y ae ^ax[2 [ (ac3Cz]>P_Q_z) & (Z2ar 6 ) & <z<(y e +y fl )/2) > ? 
where for z=:r 6 , y 6 +if C7i + Y 6 )/ 2 > 

ac:i[z] - 2 (c[l] [k] r (k,l) , A^z) ) J 

30 and r (y 6 +lrB+l)f •--» > *J > " 

<7> After tlie octagon is obtained, the curvatures of points 
in the background are set to 255 for distinguishing them from 
the points i:i clear area. 

35 3.4 Locating singularities And Analyzing Trends 

How referring to FIGS. 7 , there are three types of core 
in a fingerprint, according to the structure of the ridges 
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around them, named 'o 1 (FIG- 7a) , 'n 1 (FIG. 7b) r and "u 1 
(FIG* 7c), respectively. An 'o' core may appear in a whorl p 
an 'n' may appear in a whorl, a double loop or a loop, while 
a , u 1 may appear in a whorl, a double loop or a nodding loop, 
5 However, an/ core or delta is a point such that the 

directions of ridges around it are very inconsistent, so its 
curvature is very high and may be higher than the curvatures 
of its neighboring points* Visually, there always are many 
bright poinbs on the curvature array of a fingerprint; some 
10 of them indicate the position of a core or delta, while 
others of tiem indicate a scar, fold or noise. 

According to the publications of Sheji in 1986 P a point 
oh a digit a L image is called a singularity if it has a 
maximum curvature value that is greater than a threshold 
15 among its 8 neighboring points which form the corners and 
line midpoints of a square in which the point is centered. 

A poini P=(k,l) is called a singularity if its curvature 
is not less than its 8 neighboring points 1 and not less than 
a predetermined threshold p_ci, i.e. • 
20 1. (<?(P) > c(P+Q,)) & (c(P) £ c(P-Q t )), (1^1,2,3,4); 

2. C[P) 2t p_cl; 

There will be some singularities appearing in the region 
near a core or delta. But usually there will be some 
singularities appearing in the region near a scar, fold or 

25 noise, too. For the purpose of recognition of a core or a 

delta among singularities, analyzing the structure around a 
singularity is necessary. 

Now referring to FIGS. 8 r the ridge -flow is different 
around a core, a delta or an ordinary point. The distinction 

30 can be described easily by using the concept of ridge trend. 
A ridge trend of a point is defined as the direction of 
ridges that are near the point and run off from the point. 
There are tliree ridge trends for a delta, two trends for an 
'a 1 core or ordinary point, one trend for a l n l or ■u 1 core, 

35 and no ridge i trend for an 'o 1 core. The number of ridge 
trends of a point is called the forkedness of the point. 
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To find the ridge trends of a point P, a series of 
digital circJ.es {0 S > with P B as the center and various radii 
s are used. I'or every point P B on a circle O,, the difference , 
between the local direction at P s and the direction of vector 
5 PpP^ is calculated and stored in an array dd[]. The 

difference when P C P 3 extends substantially in the direction 
of a ridge trend will be very small, as shown by the curve 
minima FIGS.S-, so the ridge trends will be decided by finding 
all ™««i«nn« in dd[], FIG- 9a shows the pattern of array dd[] 

10 around an 'o' core, FIG. 9b around an 'n' core, FIG. 9C 
around an 'a' core, FIG. 9d around a delta. Fig. 8e 
represents a noise point or a scar. Fig. 9e shows the 
pattern of dcl[j around the point in Fig. ae. 

Fourier transform and reverse transform are used on dd[] 

15 for reducing the effect of noise. The forfcedness can be 

decided by tie power spectrum,, while the trends can be found 
on filtered c\d[]. 

The metfcod for determining the forkedness of a point 
will be described with reference to the power spectrum of a 

20 Fourier transform FT: o[0J, where j is the order 

of Fourier tl.e transform. Firstly, if o[0] is the maximum of 
all otjl an* 0 t°] is not smaller than a predetermined 
threshold p_> 0 , then it is an «o« core. Secondly, if «[13 is 
the maximum cf all otj] except e[0] f and «[l]2P_W 1r then it 

25 is a 'n' or 'u' core. Else if e£2] is the maximum of all 

except e[0] and «[2]&p_o 2 , then it is an »a' core. Else 
if o[3] is tte maximum except tt[0] and <a[3]fep_«o 3 , then it is 
a delta point:. Here p_o 0 / P-«V P-^z 5Jid P- 6 ** are a11 
predetexminec. parameters. 

30 m the case of Fig. 8e, o[2] to be the maximum, but it 

is not an »a» core. TO determine whether it is an »a' core 
is depended en its trends. If the difference of the two 
trends is veiy large, for example, it is great than*2r/3, 
then it would, not be an 'a» core and would be ignored. 

35 However it is notable that both the trend directions and 

forfcedness at. a point are due to the radius of the analysis 
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circle. The algorithm for analyzing the trends of a point P 
is as follows: 

<1> L6tS = S nin , 

Where is the minimum radius of digital circles for 
5 trend analysis and is the maximum* 
<2> Calculate dd[], 

if (s > s^) then reject and return; 
else {for every point P 0f on digital circle 0 & do: 
<dd[i] = |d(P 8f ) - dv(P 6< -P)| % 7r; 
10 if (dd[i] > ir/2 then dd[i] - sr-dd[i] ; } 

ns «= #{X | (x on O s ) & <c(X) > p_c 2 )}; 
if (nz > p_nj) , then (s = s+lj goto <2>} 
else goto <3>;} 
i.e. if the number of points, p fl ., on o a with. c(P B1 )>p_C2 is 
15 more than p.ji 2 , then increase the radii s and recalculate 
dd[]. 

<3> Derive Fourier trains form on array dd[]* 

a[j] - S (dd[i] *cos(2-?r-i' j/m) , i, (0, . . • ,m-l})/m; 
b[j] ~S (dd[i] •sin(2-*-i- j/m) r i, { a, . . . ,m-l))/m; 
20 for j = 0, • . . , n; 

where m=#0^, 

The poirer spectrum is 
s=s win -l; - (atj])*+(*[j]) 2 ; for j-o, n. 

where n is i±e order of the Fourier transform. 
25 if (o[0] > p_w 0 ) (P is an l o» core, return;] 

else {Aet k satisfying, 

«;k] -maac{«[l] r a[2], o[3]}? 
i:: (u[kj < p_» k ) then goto <2>; . 
el.se if k is equal to 1 then 
30. (calculating the trend direction, 

td = arctan(b[l]/a[l]) <-tt < td < *r) ; 
if (td > 0) then P is an 'n' core; 
else P is a f u f core; 
return ? 
35 } 

else make the reverse Fourier transform; 
(dd[i] = £(a[j] •cos(2-?r-i> j/m)- 
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10 



b[j] •sin(2-iT'i'j/Bi) , j r (0,i,...,n)) ; 
} 

} 

where p_.Oj(j-O r . . . ,3) are predetermined parameters* 
<4> Finding each minimum value ddcij] in dd[][, i« e - for 
each dd[ij] tbat satisfies 

ddCijlSdacCij-l)^] * dd[i j ]£dd[(i i +l)%a & dd[ij]<p_dd, 

for j=i,« . -*1« 
if (1 is; not equal to k) then goto <2>; 
else {tie forkedness at P is k; 

thei trends at P are dv(P fiTJ -P) , j=l, ...,k- 

) 

return; 

where p__dd is a predetermined threshold for finding minimum 

IS points. ddC(i r l)%m] and dd[ (ij+l) %m} are the values of dd[] 
at neighborirjg points of P 8ji on the circle. 

In the equations presented above, n is a constant that 
is much smaller than m, the parameters s Q , n and are 
predetermined.; P _ C , P _n 2 and a. (j=0,...3) are thresholds. 

20 More specifically, p.r^ is selected to eliminate curvature 
values which axe so large in number that the associated 
direction value is unreliable. The threshold p_c is selected 
to obtain points with a suitably high curvature value. 

According to preferred embodiments of the invention, s ofn 

25 may be equal to 5 and may be equal to 20. 

The procedure for finding all cores and deltas of a 
fingerprint €*oept the 'a» core is as follows: 
<1> 2111 singularities i.e. maximums of curvature array in 
the clear region, are located j 

30 <2> Every singularity is analyzed for finding trends by the 
above algorithm. 

<3> select cne point in every set composed of the same kind 
of cores or Celtas located together as the representative. 
•axe. criterior. for selecting -o' core is »[0], for -n' or »u« 
35 core is e[l], for delta is o[3]. 

in section 3.7, a method will be provided for finding 

the 'a' core. 
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3.5 Line tracing 

A digital curve in a fingerprint is called a contour 
line if it is in keeping with the local ridge direction at 
every point on it. Usually, a contour line can be obtained by 
5 starting from a point on the fingerprint with a trend of the 
point as the initial direction and extending or tracing 
progressively according to the direction array. Especially, 
if the initial point is a core or delta of the fingerprint, 
then the extended trace is called a shape line of the 
io fingerprint _ 

In the method provided below for accurately tracing 
contour lines, every tracing step moves by just one point; 
meanwhile ttie errors of each coordinate and direction are all 
accumulated for correcting the tracing; furthermore the 
15 tracing will stop at the right place. Where JcO and 10 are 

the column and row with minimum curvature in the clear region 
produced as described previously, cl, ck and cd are the 
current row r column and direction values, the algorithm for 
tracing a contour line starting from point (£0,10) with 
20 initial direction do is as follows: 

dl r dfc» dd are the accumulating differences of 
coordinates and direction in tracing respectively; tll[], 
tlk[], tld[ | are the arrays of coordinate and direction of 
referring points in tracing, di is a parameter , ac is the 
25 average curvatures in the segment of a line. 
<1> Initializing the variables, 
dl=d]o»ad=fcO; cfc=k0; cd=dO; 
1=1; a<:=0. 

tll[0]=;lO; tlk[0]=K07 td[O]=d0; 
30 goto<2:'. 

<2> Stepping to next point , 

while icd-tld[i-l] < -jt) cd « cd+ir; 
while icd-tld[i-l] £ sr) cd = cd-ar; 
if (|ccl-d0| > p_d3 then goto <3>, 
35 where p_d 3 is the limitation of the difference between 
current direction cd and initial direction do. 
ac = ac+c[cl] [ck] ; 
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if (i > p_l) 

{ac=ac-c[tll[i-di] [tUc[i-di] ] ; 
if (ac > p_ac-p_l) then goto <3>; 
where p_ac is a threshold. 
5 dd = cd-tld[i-l] ; 

if (|dd| > p_d2) then goto <3>? 
where p_d2 is the limitation of accumulative difference of 

direction. 

if (dd > p_dl) then 
10 [dd > p_dl; 

cd = cd— dd 

> 

else if (dd < -p_dl) then 
{dd = dd+p_dl? 
3_5 cd cd— dd; 

} 

else dd = 0; 

where pjdl means the maximum value for correcting the 
direction. 2 he increase of cl and ck depend on sin(cd) and 

20 cos(cd): 

if (|sii;Ccd)| < lcos(cd)|> then 
tcJ - cl+sign(cae(cd) ) f 
dk = dk^tan(cd) ; 
if (|dx| a i) then 

25 

{ck ~ ck+sign(dk) f 
dk = dk-sign(dk) , 

) 

} 

30 else then 

(c>: = ck+sign(sin(cd)) r 

dl = dl+ctan(cd); 

if (|dl| £ 1) then 

{cl = el+sign(dl) , 
35 dl - dl-sign(dl) , 

} 

} 
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if point (ck, cl) is out of the clear region, then goto 

<3>; 

else, saving the current coordinates and direction, 
<t.U[i]=cl, tlk[i]=c*, tld[i]=cd, i=»i+l; 
5 cci » d[cl][ck]; 

goto <2>i 
} 

<3> Determining the length of the traced line, 
i = i-1, 

10 if <e[t,ll[i]][tlk[i]] < p_ac) then goto <4>, 

else (i = i-1? 

ij: (i > 0) repeat <3>, v 

else goto <4>; 

J 

15 <4> the teiced line is {tll[j], tlk[j], tld[j], j=t),i, - ■ > ,i; 
xciturn, 

where di, p..ac r p.d^ p_dg and p_<% are all predetermined 
parameters. 

A contour tracing will be stopped if one or more of following 
20 conditions j.s true: 

(D current point is out of the clear region. 

(2) Avorage curvature of last several points is too 

high. 

(3) Rol:ated angle from initial direction is too large. 
25 A simil.ar algorithm is used to trace any normal line of 

a fingerprint; the normal line is defined as a digital curve 
on the imag<* that is perpendicular to the local ridge 
direction everywhere. The algorithm can be obtained from the 
one above by replacing the assignment cd=d[cl][ck] in step 
30 <2> With cd ! =d[cl] [ck]+7r/2 or cd=d[clj [ck3-7r/2., The normal 
lines are uised in a novel method described in section 3.6, 
below, for inacroscopically locating the center of a 
fingerprint. 

Line tracing is a basic algorithm in this fingerprint 
35 processing jsystem. It is important in locating the 

coordinate iaxes, extracting the shape features and detecting 
minutiae etc. that will be described in following sections. 
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The linci tracing operation described in this Section is 
used to trace* shape lines from the center of a delta. 

3.6 Macrosccipic Method For Locating Coordinate Axes 
5 The typsi or shape of a fingerprint, especially if it is 

characterized' by a small whorl or loop, may be ambiguous due 
to distortion and noise produced when the impression is 
taken. Some depressions of a small whorl look like a loop or 
tent arch, wtiile some impressions of a small loop look like a 

10 plain arch- Therefore, there must he some consistency in the 
rules for deciding the coordinate axes, i.e* the center and 
central orientation, of various types of fingerprints . 

In section 3.4, there was described a method for 
locating cores and deltas, except the "a 1 core,, of 

15 fingerprints by analyzing all singularities of the 

fingerprint. However, to find an r a' core and to determine 
the center and central orientation of a plain arch that are 
consistent with other similar fingerprints, a method which 
involves analyzing the macroscopic structure of fingerprint 

20 ridges is needed. 

Referring to PIGS. 10 again, many fingerprint types or 
shapes, for example whorls, double loops, loops, plain arcdl, 
tent arch etc*, are shown. In fact, the principal 
distinctions among them are always at the central parts of 

25 the fingerprints, while the peripheries of fingerprints are 
all very similar. Generally, at the central area of any 
fingerprint, the ridges at the upper part will form a vault, 
the ridges at left and right sides will run off from the 
central part, and the ridges at the lower part will be always 

30 plain, as shown In FIG. 11. 

Referring to FIG. 12, any core of a fingerprint, except a 
8 u' core, is always at the most curved region of the ridges 
below the vault formed by upper ridges. Generally for any 
fingerprint except a nodding loop, the normal lines will, 

35 starting from the upper part and going down, all concentrate 
together at a central region where there are most curved 
ridges of the fingerprint, i.e. where there is an 'o' core or 
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an 'n* core for a whorl or loop, or an 'a 1 core for a plain 
arch. Some of the normal lines may end at the central region, 
while other* i may obviously curve at the central region, so 
the core can be located by analyzing the singularities near 
5 each end poijit and the most curved point on' the normal lined* 
This method is very important for a plain arch, because it is 
usually difficult to determine, or locate, the center of a 
plain arch. The algorithm for locating microscopically any 
core except a 'u' core of a fingerprint is: 

10 Let 11 be the upper line border, 12 the lower line 

border, kl 1:be left column border and lc2 the r±ght border, 
<1> initializing, let the initial start point (ko, 10) for 
tracing be 

10 « 13,; kO - (k2-kl)/2; 

15 <2> selecting current start point (kO, 10), 
10 » l(i+l; 

if (c[l0][k0J > p_c3 then goto <2>; 
dk = p„k'sin(d[10][k0]); 

if (|d>:| > p_d4) then {kO = kO-dk; goto <2>;} 
20 else, t;he start point is (kO, 10) r 

do = d|i0][k0]; 
goto <:<>. 
where p^o3 <ind p_k are parameters, 
<3> Finding a vault line. 
25 A vault line can be considered as a combination of two 

contour lincts, i„e. a right contour line and a left contour 
line that both start from the middle of the print. So 
firstly the two contour lines should be traced. Starting 
from (10, kO) , two contour lines can he obtained .by tracing 
30 . with directions dO and d0+5T respectively. These two contour 
lines are then combined into a vault line* 

If the vault is not perfect, i.e. if its length is too 
short or its; chord is too slanting, or the curvature of the 
vault is too high, then goto <2>, 
35 . . else U± <vl[j] r vk[j], vd[j], j^0,.,-,lv} be the vault 
.line; where vl[j], vk[j], vdtj] are the y coordinate, the x 
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coordinate and the local direction, respectively, at the jtfa. 
point on the line and lv is the length of . the vault line. 

let i=0; goto <4>? 
<4> Tracing a normal line Joy the previous algorithm with 
5 starting point (Ic0,lO) and direction dO that 
jco = vk[i]f 10 - vl[ij; do = vd[i]-Hr/2; 
if (i < lv) then 
(i = i+p_g; 
goto <&>, 

10 ) 

else goto <5>; 

where p_d4 is a threshold for limiting the local direction of 
starting point, p_g is the gap between two normal lines at 
starting points. 

15 <5> Determine the areas of concentration of the normal lines 
(see PIG. 12) , then analy2e the singularities in the area to 
find the 'o', 'n' and 'a' core or others by the forKedness 
with rite algorithm described in 3.4. 

In the above algorithm, p_c is a threshold, p_k is a 

20 constant, in the trend analyzing of above points, if the 

fortasdness is 0 or 1, then the point is an 'o- or -n' core, 
while if the number is 2, then it is an 'a' core. 

By singularity analysis, there may be more than one 'a' 
core in the central region of the fingerprint. The criterion 

25 for selecting the most representative one among them is the 
angle difference dd between the two trends (di and da) of a 

singularity, i.e. 

dd = min( | d2-dl | , 2 • jt- | d2-di | > . 

For example, if dl=*/4 and 02=3*74, then dd=ir/2> if dl=2ff 
30 and d2=3?r/4 f then dd«3»r/4. 

The 'a' core which has the smallest angle difference 
will be selected as the most representative one. 

Referring to FIG. 13, there are two trends of an -a' 
core, one is towards left while the other is towards the 
35 right. The main trend of an -a- core is defined as the trend 
at the core side where the gap between two adjacent contour 



PAGE 512 1 RCVD AT 6/29/2005 1 :40:37 PM [Eastern Daylight Tone] ' SVR:USPT0-EFXRF-1/3 ' DNIS:8729306 ' CSID:631 549 0404 * DURATION (mm-$$):23-48 



06/29/2005 13:56 FAX 631 549 0404 



STRIKER & STRIKER 



-* US PTO 



0)054 



WO 93/07584 PCT/U592/0S446 

35 

lines is wiler than at the other side. By this rule, the 
central orientation of a plain arch is consistent with loops. 

After ihe trend analyzing for each singularity and 
locating th-a 'a 1 core macroscopically, the center and central 
* 5 orientation of a fingerprint can be decided as fallows 

sequentially: 

<1> If the; re is an l o l core in the central region, then the 
pattern mus" toe a whorl, the position of the center is the 
center of t3ie core, the central orientation is i.e. down 

10 forward, and 0° is horizontal to the right. 

<2> If, in the central region, there is an 'n 1 core and a 
»u' core, the pattern is a whorl; if there is an f n' core, no 
'u T core and more than one delta, the pattern is a whorl f and 
if there is an 'n 1 core, no 'u 1 core and not more than one 

15 delta, the pattern is a loop- The position of the center is 
the same as the center of the core r and the central 
orientation is the trend of the core. 

<3> Else i:; there is a i u t core, then the pattern is a 
nodding loop, the position of the center is the same as the 
20 center of the core, anS the central orientation is the trend 
of the core, 

<4> Else i:! there is an 'a' core in the central region, then 
the pattern is a plain arch, the position of the center is 
the same as the center of the core, and the central 

25 orientation is the main trend of the core. 

The center and central orientation of a whorl, a loop 
and an arch,, decided by the macroscopic method are shown in 
FIGS » 19, 20 and 21, respectively. These FIGS- depict tracing 
lines which have been generated to be perpendicular to the 

3D local ridge directions in the vicinity of the center of the 
fingerprint pattern - 

3.7 Shape Features And Classification 

All of the shape lines of a fingerprint can be obtained 
35 by the abov<* algorithm starting from the center of delta of 
the fingerprint with each trend as the initial direction. 
Line iiracincr is performed as described in section 3*5, above- 
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Various shape lines of fingerprints, much as shown in 
FIGS. 10, can serve to describe, the shapes of fingerprints 
accurately. According to the structural relations of the 
shape lines, the fingerprints can be classified into 18 types 
5 each with a respective topological structure. A finer 
classification may be based on the shape features defined 
below. 

FIG- 14 illustrates a technique for extracting the shape 
features of a left loop. Where c is the center, P„ is the 
10 delta, sl T , sj^ and si, are shape lines starting from P ft . 

P,.-.P 7 are points selected on sl 2 or sl 3 . The algorithm for 
extracting the shape features of a loop is as follows: 
<1> Determine center C-(3cc,lc) by the algorithm in Section 
3.6? 

3j> < 2 > Determine delta centered at P„ -(k„,l») by the algorithm 
in Section 3. i; 

<3> Trace these shape lines sl v sl 2 and si, starting from Po 
with three tnands of the delta as initial directions 
respectively - 

20 <4> selectiru j seven points P,»Ck,-,l f )» (i=l, - . . ,7) on sl 2 and 
slj, such that. 



where d o =dv(P 0 -C) . Because P 6 is always very far from C or 
beyond the bo:zder of the image, no feature will be defined by 

25 referring to 

<5> calculate the distances between c and P f , i.e. 



<6> Ridge counting between C and P 0 , 

Let £CWo>' OV Y n )} be the straight line from c 

so to P e , where i^y^frcAo) , •CWn>-f k .' 1 o) > let 
g[i| «= fLy,!^], (i=0,...,n); 
then the ridg^a counting is defined as, 
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rc(C,p 8 ) - #(g[i] I g[i-i]>g[i]<g[i+iJ» (i=i/.-n-i)). 
<7> Total cf IB shape features are defined, they are 7 
distances |*j-CJ, one ridge count rc(C,P e ), and 10 direction 
values referring to d, including the central orientation, 3 
5 trends of P„ and 6 local directions d(P { ) , (i^l, . . .5,7) - 

The 18 shape features of a right loop are extracted in a 
manner similar to the left loop. A whorl or double loop can 
be considered as composed of two loops, i.e. one left loop 
formed by tt.e left delta and center as well as one right loop 

10 formed by ti.e right delta and center, so it has both 18 left 
loop shape features and 18 right loop shape features, i.e. 36 
shape features in all. For a tent arch, the features 
referring to point P 7 are not extracted because it may appear 
at either leift or right. However, for the purpose of 

15 consistency between whorls and other shapes, a total of 36 
shape features are supposed for any fingerprint. If some of 
the 36 shapsi features can not be obtained due to noise, 
imperfections in, or the shape of, the fingerprint, these 
features ar^ each assigned a value of -1. There is no 

20 meaningful ishape feature for a plain arch, in other words all 
shape features of a plain arch are equal to -1. 

According to the structural relations of position and 
surrounding etc. among shape lines, fingerprints are 
classified :.nto 18 classes with respectively different 

25 topological structure as shown in FIGS . 10 which show 11 

whorls, 4 loops, one accident, one tent arch and one plain 
arch. Every class of whorl, loop and tent arch can be further 
classified according to the shape features. 

30 3.8 Global features and Global difference.'' 

Tne shiipe features for describing the pattern of a 
fingerprint are all based on both the center and delta.. So 
they nay be affected by the imperfections of fingerprint or 
by the nois* which distorts the center, delta or shape lines. 

35 m particulitr, there is no shape feature defined for plain 

arch, so thit for the purpose of practicality and consistency 
of fingerprLnt system, the features for describing the 
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pattern of a partial, noisy or plain arch fingerprint should 

be considered. 

One of liie most important parts in this invention is a 
method for defining and extracting' the global features of 
5 various fingisrprints to represent their pattern naturally and 
consistently by referring to the local ridge directions. 

Generally , the global features of a fingerprint provide 
a basis method for representing the ridge direction array of 
the fingerpr;jit. These features must be obtainable for any 
10 kind of fingerprint, and be effective in pattern matching of 
fingerprints. 

The simplest method for defining global features is to 
select some points on the direction array and take the local 
ridge direct:.ons at each point as features. So if the amount 
15 of points is large enough, then the accuracy of repression 
will be fine enough. Especially, as shown in KCGS. 15, the 
points can b<» selected to form a circular, or polar, array or 
a rectangular array. 

First method, the points can be selected on several 
circles with a common center- Referring to Pig. 15A, 
c-CK^l«> is H =a a center of a fingerprint, d, is the central 
orientation, or direction. There are n circles 0,- 
(i-0,. ,n-lj # vith a common center C and different radii r, 
(i=o,...,n-l] . There are m t selected points F f j-(k,j,i t j) 
25 (j=o, ...,1^-1) on o f segmenting the circle equally 

(i=*,...,n-i) . The global features gf of a fingerprint are 
defined as a set: 

gf==[gff ((kutlfj)) | j=o,...,mi-i? i=i,. ..,n-i]; 

where 

30 k,j:=int(k 0 +r f • cos (d^+j • 2 - jr/m,) +0.5) ; 

l t fsint (l0+r { • sin Cd 0 +j ■ 2 • ir/m, ) +0 . 5) ; 
255, if c(P) > P_C4> 
gffC*) - < int[d(P) ■p_jr/»r+0.5)%p_ir, elsewhere. 
35 i.e., gfCP) equal to 235 when the ridge directions around 
P are not cl«!ar or the curvature c(P) is greater than a 
predetermined threshold p_c4; else where gff (F) represents 
the direction d(P) in one byte with a value between 0 and 
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251; parameter p_ir is selected to transform the range of 
angular value from [0,ir] to [0,p_wj for storing it in one 
byte and reserving enough accuracy. 

In an «onbodiment, the values of parameters are n~9 and 

5 3o f =64 (i=a, ,n), so there are total 9- 64-57 6 points 

selected, at id the global features of a fingerprint are 
composed by 576 bytes. If the number of points selected on 
each circle o, is equal to others, then gf can be simply 
stored in an array of bytes: 

gf=[ff£[£] If} (=ffff[P S j) ) J:7=0, . . .,n?-l;i=0, . • . ,:n-l;} 

10 

For exsjaple, the global features of an arch, loop and 
whorl are shown in FIGS. 22, 23 and 24, respect ively , where 
the center of each pattern is at the common center of the _ 
concentric circles and the central orientation of each 

IB pattern is represented by a short line extending from the 
center. Th^re are 9 circles in each image, and each circle 
is composed of 32 points (for showing more clear than 64 
points) . Fc.r each point selected, if it is not in background 
and its curvature is not high, then the local direction is 

20 represented by a line centered on the selected point; the 
remaining selected points are each represented by a dot r as 
is particularly apparent at the bottom and lower portions of 
the left-hand and right-hand edges of fig. 22. 

Ofte difference between two shapes of fingerprints always 

25 will reflect on their direction arrays* So that it also 
would reflect on the global features which represent the 
direction arrays. For this purpose, an important 
measurement, called global difference between two sets of 
global features, is necessary. 

30 In the case of Equation (9) , the global difference gdi 

between two fingerprints by their global features gfl and gf2 
is defined c.s: 

gdi(gfi, gf2) min{S(f_dg(dl-gflfi] [j]/ 
d2-gf2ti][(j+r)%m])/#Ml(r), (i,j), Ml(r)), 

35 . r='Q,...,Bi-l* #Ml(r)>0}. 
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where dl and d2 are central orientations of two fingerprints 
respectively, set Ml(r) means 

MlCr)^{ | gfl[i][jJ<P-»" * gf2[i][(j+r)%mj<p_ff); 

function t_dg() means 
5 f_dg(x,y) = f_i(min(|x-y|%p_7T/ p_Jr-(|x-y|%p_»") ) ) ? 

where f_i(z) is an increase function of z, in the 
embodiment, f _i(z) = a** So tHe global difference between 
two fingerprints is calculated by matching their global 
features with various radii r to find the minimum difference. 
10 Second mathod, the points can be selected on a grid with 

n rows and m columns, as in FIG. 15b r so the global features 
can be stored in an array gf [] [J such that: 
gf ti] [j ]=d[Y« +i ' d yJ [fc^+j >dxj ; 

i«0, . „ . ,n-l; j«0 , . , , f m-1. 
15 where x^ and jr^ are the coordinates of the left upper corner 
point on the iirray, dy and dx are increases of row and column 
respectively . 

The global difference gd2 between two fingerprints by 
their global features gfl and gf2 is defined as: 
20 gd2(gfl r crf2)^ain{2:(f_dg(dl-gfl[i3PCj]/ 

d2-<rf2[lHk]/*M2<l/fc) r (i,i)r M2(l,k)), 
complete 1=0 , * • • , n— l, 3c=0 r * • . , m— 1- 

The globiLl difference can be used for finely classifying 
fingerprints ::ti a database, or selecting similar fingerprints 

25 in a database to reduce the difficulty of minutia matching 
during a searching procedure. 
3.9 Detecting Minutia From Gray Level Image 

Minutiae are very important traditional features of a 
fingerprint, Eind are used in final verification of identity 

30 of two fingerprints. Usually minutiae are described with 

respect to the: pattern of fingerprint ridges. There are many 
types of minuiiaa on a fingerprint, for example as shown in 
FIGS- 16, endings (a), bifurcations (b) , islands (c) r eyes 
(d), bridges |e), etc. In brief, minutiae are singularities 

35 of ridges* 

However ridges always coexist with valleys on a 
fingerprint, 5Jid each feature or minutia of ridges always 
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corresponds to a change in valleys , so that minutiae can be 
described in terms of valleys, too. In general r an ending of 
a ridge is e. bifurcation of valleys while a bifurcation of 
ridges is an ending of a valley, an island of a ridge is an 
5 eye of a valley and an eye of a ridge is an island of a 

valley. Refcxring to FIGS, 17, the exceptions may appear at 
cores and deltas. The description of minutiae that are just 
at a core or delta in terms of ridges is different from the 
description in terms of vaileys. However, the descriptions of 

10 minutiae should be consistent by being all in terms of 
valleys or all in terms of ridges. 

For automatic detection of minutiae, the novel method 
provided here is based upon tracing the valleys rather than 
the ordinary method which is based upon bin arizing, thinning 

15 and smoothing the ridges. Generally in a fingerprint image, 
the quality of valleys is much better than that of ridges, 
primarily because of the following reasons: firstly, there 
are no sweat glands in valleys; secondly, the widths of 
valleys are more even than those of ridges; and thirdly, the 

20 gray levels in valleys are more even than in ridges. Although 
there will be incipient ridges in the valleys of some 
fingerprints that may affect valley tracing, all ridges of 
every fingerprint have sweat glands that may affect ridge 
tracing. So in general, the result of valley tracing should 

25 be much better than ridge tracing. 

The algorithm for tracing a valley with an initial point 
(k0,lO) and direction dO is similar to tracing a line, except 
that it uses a key technique that keeps the step points in 
the valley. 

30 . :Let ftl][k] be an image array. Its element f [1] [K] 

equals the gray scale value of a point (k,l) , it will be set 
to ~1 after it has been traced, ag is the summation of gray 
scale values of the last p_l points in a tracing line. Array 
tlg[] is used to store gray scale values of traced points. 

35 The definitions of other variables are the same as for the 
algorithm of line tracing in section 3.5. 
<1> Initializing the variables. 
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<U=dfc=acl-o; 

cl='lO; cJc=kO; cd=dO; 
i=l? ac= r O- 

tll[0]=10; fclk[0]=ko; tld[0]=do; 
5 tlg[0]=f[iO]lJcO]; 
goto <2? • 
<2> Btep to the next point, 

Accumulate curvatures of every point in the valley, 
ac=ac+c[cl]tck]y 
10 if (i^P_.D then 

{ac=ae-c[tH[i-p_l] ] [tlk3c[i-p_l] ] ? 
if (ac > p_acp_l) then goto <3>? 

> 

Accumulate gray scale values of every point in the 
is valley. 

ag*=ag+f [cl] [ck] ; 
if (i > P_l) then 

{ ag =*ag-tlg [ i-p_l ] ? 

if fag < p_ag»p_l) then goto <3>; 

20 } 

if (cd-tld[i-lj £ -jt) then cd = cd+TT? 
if (cd-tldCi-l] £ ^) then cd = cd-?r; 
if (|cd-do| > p_d3> then goto <3>, 
where. p_d3 is the limitation for changing current direction 
25 cd per step in tracing* 
ac=ac+c[d] [ck] r 
if fi > p_l) then 

{ac=ac-c[tll£i-p_l] J [tlkf i-p_l] ] ; 
if (ac > p_ac a p_l) then goto <3>; 
30 where p_ac is a threshold. 

} 

dd*=cd-tli[i-l] ? 

if (|dd| > P_02) then goto <3>; 
vherfc p_d2 is the limitation of accumulative difference of 
35 direction. 

if (dd > p_dl) then 
{dd = dd-p_dl; 
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c i=cd-dd 
) 

else if (dd < -p_dl) then 
{dd^dd+p_di; 
5 oi=cd-dd; 

} 

else dd^O; 

where p_dl neans the maximum value for correcting the 
direction. The increase of cl and ck depend on sin(cd) and 
10 cos(cd>; 

if (|s:Ln(cd)| < |cos(cd)[) then 
{ c:l=cl+sign (cos (cd) ) , 
d]:=dk+tan)cd) ; 
ijr (|dk| S 1) then 
15 (ck=ck+sign(dk) . 

dk^dk-sign)dk) , 
} 

) 

if any point (x,y) of ((dc,ci), (kl,ll) r (kr r lr)> is out 
20 of the cleat* region or f[y][x]<0, then goto <3>? 

else save current coordinate, direction and gray level 
values; 

{tll[i]=cl; tlk[i]=ck; tld[i]=cd; i*=i+l; 
f[cl]tck]=-l; 
25 cd^dCcinck]? 

dd«dd+(f [ll]tkl]-f[lr][krj) -p_ga,- 
goto <2>; 
) • 

Where (kl,ll) and (kr,lr) shown in FIGS, IB are called the 
30 left point and right point of current point (ck,cl) 

respectively. Both of them are the 4 -neighboring points of 
current point and 8 -neighboring points of the previous point 
(tlk[i] r tll[i]> . P_ga is a predetermined parameter of 
modifying direction by difference of gray scale values, 
35 <3> ' Determine the length of traced valley , 
i=i-l; 
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if (Cc[fcll[i]] [tlk[i3]<p_ac) & tlg[i]>P_ag)> then goto 

<4>, 

else (f [tll[i] ] [tlk[i] ]=tlg[i] ? 

5 if (i>0| then repeat <3>, 

el;3& goto <4>. 
} 

<4> The traced valley is (tll[j ] , tlk[ j ] , tld[ j] , 
(j-0,1,. . .,!)], 
10 return? 

Where p._l is a constant, p_ac, p_dl and p_d2 are all 
thresholds - 

Referring to FXGS-lS r both the left point (£1,11) and 
the right po;Lnt (kr,lr) are 4-neighbqring for current point 

15 (ck,cl) ans n-neighboring for previous point 

(tlk£i-ll,ti:.[i-i]). (ck,cl) may be replaced by its 
4-neighborincr point (kl r ll) or (kl,ll) according to their 
gray level ♦ 23iis algorithm is similar to that for line 
tracing with gray level as an additional factor. 

20 In detail/ the tracing line will firstly step from the 

prior point to the current point temporally, then the gray 
scale of two 4_jieighboring points, i.e. (kl,ll) and (kr,lr), 
Of the current point are considered. A point is selected to 
be a valley point if its gray scale value is higher than or 

25 equal to the other points in the neighborhood. 

A valley tracing will be stopped if one or more of 
following conditions is true: 

(1) The current point is out of the clear region. 

(2) The average curvature ac of the last p_l points in 
30 the tracing is very high, i.e. greater than p_ac. 

(3) Any previous valley trace is touched • 

(4) The average gray level ag of the last p_l points in 
the tracing i& very low, i-e, less than p_ag. 

The algorithm for detecting minutiae from a gray level 
35 image by valley tracing is as follows: 
<1> Let gp=cap2; 
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<2> The stait point P for valley tracing should satisfy each 
of fallowing conditions: 

(1) P is a maximum point in a 3x3 neighborhood on gray 
scale image F and f (P)>p_f* 

5 A max iiium point in a gray scale image means it is one 

whose gray scale value is not less than that of each of its 8 
neighboring points ♦ 

(2) P is in the clear region of ridges and c(P)<p_C4 
Shis means the curvature at p is smaller than p_c. 

10 (3) There is no traced line at directions c(P)+jt/2 and 

c(PJ-tt/2 within distance gp. 

Por every such point P, the valley is traced in two 
initial directions c(P) and c(P)+*r, respectively. The 
minutiae are detected as a result of the conditions for 

15 stopping a trace: If the trace is stopped due to condition 
(3), then a valley bifurcation is found; While if the- trace 
is stopped due to condition (4) , then a valley ending has 
been found. 

<3> Connect any two terminals of traces (a terminal being 
20 the start or end of a trace) ift 

(1) The two last directions of the traces are opposite 
to one another; 

(2) The positions of two terminals of traces are very 
close ? 

25 (3) Th»i average gray level between them is higher than 

p_ag* 

<4> gp=gp-l ' if gp>gapl, then goto <2>; 

. where parameters p_f , t, gapl and gap2 are all 

experimentally determined constants. Each minutia found is 
30 described by three attributes, i.e. x, y and 6_ The 

coordinates and y are same as the position of a trace . 

terminal, thct direction 8 is equal to the one of d[y] [x] and 

d[y] [x]+ff wh:.ch is closer to the last direction in the trace. 
An example of extracting minutiae by tracing valleys on 
35 a gray level of a fingerprint is shown in FIG* 25, where the 

fingerprint is the the same as in FIGS* 21 and 22 and the 

gray level if; reversed. 
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3*10 Quality Level And Vector 

The features of a fingerprint may be effected by many 
factors , for example noise level, the effective area of the 
clear region, the position of the center, the number of 
5 minutiae , and so on. sometimes the factors are due to the 
quality of tee finger itself, while at other times they are 
due to the impression or input device. A quality level qJL 
should be provided after image processing in order to make 
possible an automated or operator controlled decision as to 
10 whether to accept, reject, or reinput the fingerprint image, 
or if possible to take a new fingerprint impression. The 
quality level can be described in detail by a quality vector: 

q_v=(g_n, q_a r <LPr g_ m / ; 
where each factor is calculated as follows: 
15 <1> The noise level qjl refers to the average curvature in 
the whole clear region, i.e. 

q_n ~fjg_n ( a_c) t 

where 

a_c-S(c(X), X, {X|c(X)<l})/#[X|c(X)<l}? 
20 f_q_n(z) is a a increase function for z in the range 0 to 1. 
In the embodiment q_n is defined as, 
0, when a_c<el; 
l r when a_c>cl & a_c<c2; 
q-n 2, when a^c>c2 & a_c<c3 1 

25 3, when a_c>c3- 

Where cl, C2, and c3 are all predetermined experimental 
values* If tie average curvature of an image is very small, 
then q_n-o, i.e. the image's quality is good; while if the 
average curvature is large, then the noise in the image would 
30 affect the processing, q_n will equal to 1, 2 or 3 depending 
on the noise. 

<2> The effeistive area q_a represents the number of global 
features which define a direction, in the case of Equation 

35 Q^a?-f_q_a(#{x [ (x in gf) & (*<p_Jr)})? 

where function f_q_a(z) is increasing for z between o and 
255, to ensure that q-a is 1, 2 or 3. When q_a equals 0, the 
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quality of 1:he image is good, otherwise the larger the value, 
the worse the quality. 

<3> g p depends on whether the position of the center 

(JtO,lO) is : P n the central region CR of the image, i.e. 

5 q_P=(is (kO,10) in CR)? 0 ; l r - 

Xn the embodiment, 

CII={(JC,1) | (1/4 < 1 < L-2/3) k (K/3 < 3C < K-2/3)}, 

<4> q m depends on the number and average quality of 

minutiae a_iiq, i.e* 

io o, if *nm < p_nm; 

CLm = { 

f_q_m(a_mq) , otherwise, 
where p^nm de a predetermined threshold, fjmq(z) is an 
increase furction for s« in the embodiment, p_nm=ie, 
15 f_inq(2)=z/4. 

<5> q_h means help level that presents the reliability of 
the center C by the average curvature around C, 
q_h=fjjJi(2(c(X) , X, NC)/#NC) ; 

where 

20 NC=(X J |X-C|£p_rl & |X-C|£p_r2. 

f_q_h(s) is an increase function for z. 
Finally, 

q_l=q_n + q_a + q_p + q_m + q_h; 
It is anticipated that the invention will be implemented 
25 by means of a general purpose digital computer system which 
is programmed in accordance with the algorithms described 
above and is provided with an appropriate graphics input 
device capable of scanning a fingerprint image and inputting 
gray level inage point brightness values and displaying and 
30 printing and writing the results of the image processing 
procedures ii an output device. 

In an e;ubodiment of the above method, the following 
parameters may have values in the ranges specified below: 



section Parameter Range 

35 . 3.2 r [5, 30] 

p_v [4«r*r, 60tt] 

3.3 - p_nl [1, 5} 

PJ=LP [0.5, 1] 
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3.4 


P_cl 


[0.6, 1] 




P_ M o 


[1*5, 2*9] 




P-J»i 


[0-5, 1.5] 




PJ»2 


[0,5, 1-5] 






[0-5, 1-5] 




P_C2 


[0.6 r 1] 




P_n2 


[1# io] 




p_dd 


[0»1, O.oJ 


3.3 


p_dl 


r n "i ft B 1 
[0*2/ &J 






[0*0, X.DJ 




p_d3 


— *— 

grease taon u 




P-J- 


[5, 15] 




p_ac 


[0.6, 1] 




p_C3 


[0.6# 1] 




Vjt 


[1/ 10] 




p_d4 


[10, 100] 


3.8 


p_c4 


[0.6, 1] 


3.9 


P-ag- 


[0, V] 




p_g* 


[0.01, 0.1] 


3.10 


p_nm 


[&r 20] 






[5/ 20] 




p_rl 


[1, 10] 




p_r2 


[10, 20] 


The invention thus provides 


a method for calculating the 



25 global difference between two stripe patterns toy means of 
their global features used in finer classification and 
search. 

While the description above refers to particular 
embodiments of the present invention, it will he understood 

30 that many modifications may be made without departing from, 
the spirit thereof. The accompanying claims are intended to 
cover sudh modifications as would fall within the true scope 
and spirit of the present invention. 

The presently disclosed embodiments are therefore to be 

35 considered in all respects as illustrative and not 

restrictive, the scope of the invention being indicated by 
the appended claims, rather than the foregoing description. 
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and all changes which come within the meaning and range of 
equivalency of the claims are therefore intended to be 
embraced therein* 
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WHAT IS dAHIED: 

1. an automatic method for identifying an image having 
varying brightness values constituting a pattern of stripes 
5 containing m:jitrtiae f comprising: 

dividing the image into a rectangular matrix of 
image points and providing a representation of the gray level 
value of the brightness at each image point; 

transforming the representations into a selected 

10 brightness riinge; 

calculating the average direction of local texture 
at each point on the image with a function of gradient models 
in a neighborhood of the point; 

calculating the curvature or the inconsistency of 
15 directions or. local texture at each point on the image with a 
function of crradient models in a neighborhood of the point; 

separating a useful portion of the image, which 
contains an Eiccurate representation of the pattern of 
stripes, froii noisy background; 
20 extractjjig the global features of the image to represent 

the direction array by selecting some points on the Image and 
storing the directions at these points; 

finding singularities in the useful portion by 
comparing curvature values of selected matrix points in the 

25 useful portion; 

producing representations of the shapes of selected 
stripes of the image based on the curvature and direction 
values; 

locating a standardized coordinate axis system 
30 having an origin on the image by locating an intersection of 
lines normal to selected stripes in a selected region of the 
image; 

producing representations of selected 
characteristics of the shape of one of the selected stripes 
35 on the basis of the distance between the origin of the 

coordinate aicis system and points on the selected stripe in 
selected directions from the origin and producing a 
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representation of a global feature function of the 
fingerprint based on a set of values each representing a 
difference tietween the direction between the origin of the 
coordinate and a point on a stripe and the direction of 

5 the stripe jattern at that point f 

determining the locations of minutiae in the image 
on the coordinate axis system from the transformed 
representations of the gray level values of the brightness at 
each image point; and 
10 producing a representation of the quality level of 

the image. 

2. The method defined in claim 1 Vherein the image 
represents a fingerprint. 

15 

3* The method defined in claim 2 comprising employing 
a quick recurrent algorithm for calculating the direction and 
curvature with same size as the image, by the steps of; 

(1} Calculating four gradient model arrays about 0, 45/ 
20 90 and 135 dagrees respectively; 

(2) Calculating the average gradient models in each 
' neighborhood recurrently; 

(3) Using tables instead of operations of arc- tangent, * 
wherein the curvature is a measurement of accuracy about the 

25 direction at the same point. 

4« The method defined in claim 2 comprising segmenting 
an octagonal clear region of fingerprint ridges from 
background a: id noise by means of eight straight lines 
30 according to the curvature array. 
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5. TlLe method defined in claim 2 comprising analyzing 
the ridge tarends around a point in image by means of the 
difference of local direction and vector direction at each 
point on th«t digital circles with various radius and deriving 

5 a Fourier tiransform and reverse transform for deciding the 
forkedness ly power spectrum and finding the trends by 
filtered differences. 

6. Tte method defined in claim 2 comprising performing 
10 line tracing based on direction and curvature values by 

accumulating errors of coordinates and directions to correct 
the trace fcr extracting contour lines, shape lines and 
normal lines. 

15 7- The method defined in claim 2 wherein the 

coordinate axis is located by locating the center and central 
orientation of the fingerprint macroscopically by means of a 
vault line and normal lines* 

20 8. Ths. method defined in claim 2 comprising locating 

the 'a 1 core of a plain arch and its main trend 
macroscqpicaLly that are consistent with other types of 
fingerprints . 

25 9. Tha method defined in claim 2 wherein the shape 

features for both loops and whorls are extracted from shape 
lines, the sCiape features of a whorl being composed by two 
parts of two loops referring to two deltas respectively - 

30 10* l!hii method defined in claim 2 comprising producing 

a general cl.*issif ication by means of relations among shape 
lines and fine classification by means of shape features. 
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11. Tae method defined in claim 2 comprising 
identifying global features by angle values on some points 
selecting from the direction array, the points are located 
spatially o;i a lattice or on several circles with a Common 
5 center- 



12. The method defined in claim 2 wherein the minutiae 
are located in terms of endings or bifurcations of valleys. 

10 13. The method defined in claim 12 wherein valleys are 

traced and tximitiae located on the gray level image by means 
of contour tracing referring to gray levels, comprising the 
steps of i 

(1) selecting starting point for valley tracing, 
15 (2) connecting two terminals of traces which satisfy 

certain conditions, 

(3) changing the gap of valleys. 

14.. Tl.e method defined in claim 2 wherein the quality 
20 level of the fingerprint is determined on the basis of 

quality vectors referring to the position of center, number 
of minutiae, noise level, and area of clear region in order 
to decide automatically or suggest the operator for 
acceptation, rejection or reinput of the fingerprint image. 

25 

15« The method defined in claim 2 further comprising 
calculating global difference between two fingerprints by 
means of their global features used in finer classification 
and search. 



PAGE 72/82 * RCVD AT 6/29/2005 1 :40:37 PM [Eastern Daylight Time] * SVR:USPT0-EFXRF-1/3 * DNIS:8729306 1 CSID:631 549 0404 * DURATION (mm-ss):2M8 



06/29/2005 14:01 FAX 631 549 0404 



STRIKER & STRIKER 



-+ US PTO 



@073 



WO 93/07584 PCT/U59Z/08446 

1/9 



FIG.1 



INPUT: FINGERPRINT 
GRAY LEVEL IMAGE 



PREPROCESSING 

I 



CALCUUTI) D AND C ARRAYS}- 



CLEANIMG BACKGROUND 



FINDING CORES AND DELTAS 



DECIDING COORDINATE AXIS 



I 



TRACING SHAPE LINES 

I 



EXTRACTING SHAPE AND 
GLOBAL FEATURES 



I 



DETECTING MINUTIAE 

• 1 



CHECKING QUALITY 



CLASSIFICATION 



TREND ANALYZING 



LINE TRACING \ * — I 



r OK?^ 
Y 



N 




REJECT 



OUTPUT: FEATURES OF 
FfoGERPRINT 



PAGE 1W. * RCVD AT 6/29/20D5 1 :40:37 PM [Eastern Daylight Time] 1 SVR:USPT0-EFXRF-1I3 ' DNIS:8729306 * CSID:631 549 0404 * DURATION (mm-ss):2W8 



06/29/2005 14:01 FAI 631 549 0404 STRIKER & STRIKER -♦ US PTO @074 

WO 93/07584 PCT/US92/08446 

2/9 



+ 0OO 



• • ■ ■ 

oooo 

O o_ o o 
O o. O O 

oooo 



■ 



CM 



CO 

EC 



-k • ♦ • • 

• • • • 

-k • • ■ • 

-k • • • • 



Si 



TO 
CO 



3 



-K • • • • 
-*•••» 



f CD 





EC 



CO 




::::: rg 



• • Q. • • CO 



TO 



C3 



PAGE 74/82 • RCVD AT 6/29/2005 1 :40:37 PIKI [Eastern Daylight Time] * SVR:USPT0-EFxSi/3 *^DNIsl872930G * CSID:631 549 0404 * DURATION (mm-ss):2W8 




PAGE 75/82 1 RCVD AT 6129/2005 1 :40:37 PM [Eastern Daylight fie] ' SVR:USPTO^EFXRF-1/3 * DNS729306 * CSID:631 549 0404 * DURATION (mm-ss):23^8 




PAGE 76/82 1 RCVD AT 6/29/2005 1:40:37 PM [Eastern Daylight Time] ' SVR:USPTO-EFXRF-1/3* DNIS:872930$« CSID:631 549 0404 1 DURATION (mnv$s):2W8 



06/29/2005 14:01 FAX 631 549 0404 STRIKER & STRIKER ->• US PTO II 077 

WO 93/07584 PCT/US92/08446 

5/9 




PAGE 77/82* RCVD AT 6/29/2005 1:40:37 PM [Eastern Day1ight'fim7]^sS:USP?oTFXFS 1^ * DNIS:8729306 * CSID:631 549 MM* DURATION (mm-ss): 23-48 



06/29/2005 14:02 FAX 631 549 0404 
WO 93/07584 



STRIKER & STRIKER 



-» US PTO @1 078 

PCT/USU/08446 



6/9 

FIG. 19 




FIG. 20 




Q1 SRRTITT IT= P.H=tT 

PAGE Ml 1 RCVD AT 6/29/2005 1 :40:37 PM [Eastern Daylight Time] ' SVR:USPT0-EFXRF-1/3 1 DN1S:8729306* CS1D:631 549 0404 * DURATION fynm-$s):2M8 




ci incmri rnz cucsrr . 
PAGE Ml * RCVD AT 6/2912005 1 :40:37 PM [Eastern Daylight Time] 1 SVR:USPTO-EFXRF-1/3 ' DNIS:8729306 * CSID:631 549 0404 * DURATION (mnMS):2W8 



06/29/2005 14:03 FAX 631 549 0404 



STRIKER & STRIKER 



-* US PTO 



0080 



WO 93/07584 



PCT/US92/08446 



8/9 

FIG. 23 




FIG. 24 




PAGE 80/82 * RCVD AT 6/29/2005 1 :40:37 PM [Eastern Daylight rime] * SVR:USPT0-EFXRF-1/3 * DNiS:8729306 * CSID:631 549 0404 * DURATION (mm-ss):2348 



06/29/2005 14:04 FAI 631 549 0404 STRIKER & STRIKER -» US PTO @081 

WO 93/07584 PCT/US92/0W46 



S/9 




SU8STTTCJTE 55HEET 

PAGE 81/82 ' RCVD AT 6/29/2005 1:40:37 PM [Eastern Daylight Time] * SVR:USPT0-EFXRF-1/3 ' DNIS:8729306 " CS!D:631 549 0404 * DURATION (mm-ss):2M8 



06/29/2005 14:05 FAX 631 549 0404 



STRIKER & STRIKER 



-* US PTO 



@082 



INTERN A1 ION AL SEARCH REPORT 



International Application No. 
PCTAJS92/Q3446 



A. CLASSIFICATION OF SUBJECT MATTER 
IPC(5) 

USCL J 82/5 

According to Mtenjatuinal Fa teat Classification (IPC) Pf to both national ctaasificattan and IPC 



B. FIELDS SEARCHED 



Minimum documertaiirm ecai ched (rhffrifirnrion system followed by classification symbols)' 
U.S. : 356/71; 332/4 



Documeflteticm searched othcz thin minimum documcntalkm to the extern that such documents ut included id the fields searched 



Electronic data base con&ultec during the international search (name of data base and, where practicable, search terms used) 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Category* 



A.P 



Citation of dMUment, with indication. Where appropriate, of die relevant p usages 



US,A, 4,0£.,03S (RIGANATI ET AL) 
04 APRIL 1978 
See entiie document 

US, A. 5,06 " , 162 (DRISCOLL, jr ET AL) 
19 NOVEMBER 1991 
See entiie document 

US, A, 5,14C,642 (HSU ET AL) 
18 AUGUST 1992 
See entire decu merit. 



Relevant to claim No. 



1-15 



1-15 



1-5 



~| Further document Or6 listed in the continuation of Bo* C. | | SCO patent family annex. 



Special 



of dud 



IT 



■cr 



440imc&ldcfiaiag *fce gnu nd Ab *f lbs ftfl wWcb b not cowdera] 
|o fee pvt of putumlar nit mice 

dnrmwni whkh merj ijiraif doubts go priority ttoiox» at which it 
cUed to rttaMi* tha pubtailioo dale «f O&OVWr ctatin or other 
•pccM mm ■pedficd,! 

referfaf ID tttl <Hll di^nann. tnc, cxhibiljoa «f Oltef 



rtftrtlrtfatpubUihqd p*|artt the trnnodanal filial ductoi fater than 
xhc priority* 



Liierdixviimi pa trtifad after <hc haeiudo&il fllbf dnt or pHorfry 

t*-' ' tt"1 m frranrn rhh ihr iprttffltfftfi fr^a flad \n imdfrtrrH 'V- 

piBso^ilfl or ABoty nbdcr^ftag the Iwcudot 

if^rC^iifttJtl of ptftimfsT rctevBODo; Iha tj fi i| Ml 



docuBenl vf fMitEiifatf' MCavifteci ifac- rttiHrtt IaycoISoq- pmwoi bje 

ridaai' to iuVOlvC m InVdUT/C tftcp WQC& tBC ijocuioenl b 
bnrd *itll «fie Of BOf* Q&O *peH ri p nn D mh . nu£ cxnabtoalxn 
» ft pci f o naWi l n'H^ ul 



Data of the acLual completion cf the international search 
11 DECEMBER 1992 



Name and mailing address of He \SAJh\ 
ContauoBioBsrol Fatcartaond Tnn4*eiiitefl 
Box PCT 

Wiishmgloo, D.C. 20231 
Fdcdirilo No. NOT AFPUC VBLE 



Dale of mailing of the- international search report 



Authorized officer 



45 



Telephone No- f703) 305.4762, 



[fa LEO H. BOUDREAU ^g^^^^ 



Form PCX/1SA/210 (aeeond shi«0CJwty 1992)* 

PAGE 82/82 ' RCVD AT 6729/2005 1 :40:37 PM [Eastern Daylight Time] " SVR:USPT0-EFXRF-1/3 " DNIS:8729306 * CSID:631 549 0404 * DURATION fmm-ss):2W8 



This Page is Inserted by IFW Indexing and Scanning 
Operations and is not part of the Official Record 

BEST AVAILABLE IMAGES 

• . Defective images within this document are accurate representations of the original 
documents submitted by the applicant. 

Defects in the images include but are not limited to the items checked: 

□ BLACK BORDERS 

□ IMAGE CUT OFF AT TOP, BOTTOM OR SIDES 

□ FADED TEXT OR DRAWING 

□ BLURRED OR ILLEGIBLE TEXT OR DRAWING 

□ SKEWED/SLANTED IMAGES 

^TCOLOR OR BLACK AND WHITE PHOTOGRAPHS 

□ GRAY SCALE DOCUMENTS 

&LINES OR MARKS ON ORIGINAL DOCUMENT 

□ REFERENCE(S) OR EXHIBIT(S) SUBMITTED ARE POOR QUALITY 

□ OTHER: 

IMAGES ARE BEST AVAILABLE COPY. 
As rescanning these documents will not correct the image 
problems checked, please do not report these problems to 
the IFW Image Problem Mailbox. 



